NASA  Contractor  Report  194952 


AD'-A284  369 

■iinuni 


ICASE  Report  No.  94^il 


ICASE 


ON  THE  GIBBS  PHENOMENON  V:  RECOVERING 
EXPONENTIAL  ACCURACY  FROM  COLLOCATION 
POINT  VALUES  OF  A  PIECEWISE  ANALYTIC 
FUNCTION 


David  Gottlieb 
Chi- Wang  Shu 


am  lafli 


MELECTE 

U 


Contract  NASl-19480 
July  1994 


94  9  13  6 


Institute  for  Computer  Applications  in  Science  and  Engineering 
NASA  Langley  Research  Center 
Hampton,  VA  23681-0001 

XJTIC  QUALITY  IIISPSGTSD  3 


Operated  by  Universities  Space  Research  Association 


ON  THE  GIBBS  PHENOMENON  V: 


RECOVERING  EXPONENTIAL  ACCURACY  FROM 
COLLOCATION  POINT  VALUES  OF  A  PIECEWISE  ANALYTIC  FUNCTION  » 

David  Gottlieb  and  Chi-Wang  Shu 
Division  of  Applied  Mathematics 
Brown  University 
Providence,  RI  02912 

ABSTRACT 

The  paper  presents  a  method  to  recover  exponential  accuracy  at  all  points  (including  at  the 
discontinuities  themselves),  from  the  knowledge  of  an  approximation  to  the  interpolation  polyno¬ 
mial  (or  trigonometrical  polynomial).  We  show  that  if  we  are  given  the  collocation  point  values  (or 
an  highly  accurate  approximation)  at  the  Gauss  or  Gauss-Lobatto  points,  we  can  reconstruct  an 
uniform  exponentially  convergent  approximation  to  the  function  f{x)  in  any  sub-interval  of  analyt- 
icity.  The  proof  covers  the  cases  of  Fourier,  Chebyshev,  Legendre,  and  more  general  Gegenbauer 
coUocation  methods. 


‘  Research  supported  by  AFOSR  grant  93-0090,  ARO  grant  DAAL03-91-G-0123  and  DAAH04-94-G-0205,  DARPA 
grant  N00014-91-J-4016,  NSF  grant  DMS-9211820,  NASA  grant  NAGl-1145,  and  NASA  contract  NASl-19480  while 
the  authors  were  in  residence  at  the  Institute  for  Computer  Applications  in  Science  and  Engineering  (ICASE),  NASA 
Langley  Research  Center,  Hampton,  VA  23681-0001. 


1 


1  Introduction 


We  continue  our  investigation  ol  overcoming  the  Gibbs  phenomenon,  i.e.,  recovering  pointwise 
exponential  accuracy  at  all  points  including  at  the  discontinuities  themselves,  from  the  knowledge 
of  a  spectral  partial  sum  of  a  discontinuous  but  piecewise  analytic  function.  In  previous  papers  [3], 
[4],  [5]  and  [6]  we  have  considered  the  cases  of  Galerkin  approximation,  which  means  that,  given 
an  orthogonal  basis  {</)fc(a;)}  over  [—1,1]  with  a  weight  function  in(i),  we  assume  that  the  first  N 
Galerkin  coefficients 

1 

fk  =  —  w{x)f{x)(l>k{x)dx,  \  <  k  <  N  (1.1) 

Ck 

of  a  discontinuous  but  piecewise  analytic  function  /(x),  are  known.  Here  Ck  are  normalization 
constants.  A  procedure  is  given  in  our  previous  paper.s  to  recover  uniform  exponential  accuracy, 
in  any  sub-interval  of  analyticity  of  f(x].  The  Fourier  Galerkin  case  (d>jt(x)  =  w{x)  =  1) 

with  a  single  discontinuity  is  analyzed  in  [3];  wave  resolution  issues  are  discussed  in  [4];  sub- 
inturval  recon.struction  suitable  for  multiple  discontinuities  are  analyzed  in  [5]  for  the  Fourier  case 
and  the  Legendre  case  ((i>k{x)  =  Lk{x),  tt)(x)  =  1),  and  in  [6]  for  the  general  Gegenbauer  case 
(4(x)  =  Cl!(x),  w(x}  -  (I  -  x^)f^-3). 

The  reconstruction  procedure  analyzed  in  these  papers  consists  of  two  steps: 


1.  Using  the  given  spectral  partial  sum  of  the  first  /V  terms,  to  recover  the  first  m  ~  A 
(lOgonbauor  expansion  coefllcicntfi,  ha.-d  on  a  sub-interval  [a,  6]  C  [-"1,1]  in  which  the 
function  is  assumed  analytic,  with  exponential  accuracy.  This  can  be  achieved  for  any 
L\  function,  as  long  as  we  choose  A  in  the  weight  function  of  Gegenbauer  polynomials 

to  be  proportional  to  N .  The  error  incurred  at  this  stage  is  called  the  truncation  error.  __ 

or 

2.  For  an  analytic  function  in  [a,  6],  proving  the  exponential  convergence  of  its  Gegenbauer 
expansion,  when  the  paranmter  A  in  the  weight  function  is  proportional  to  the  number 

m. 

of  terms  retained  in  the  expansion.  The  error  at  this  stage  is  labeled  the  regularization  - 


!  Availability  Codai 


error. 


LflUU 


In  this  paper,  we  will  consider  the  case  of  collocation.  We  assume  that  the  point  values 
where  a:,  are  the  Gauss  or  Gauss-Lobatto  points  of  the  orthogonal  basis  {<^fc(x)},  of  a  discontinuous 
but  piecewise  analytic  function  /(x),  are  given.  Alternatively,  we  assume  that  we  have  an  accurate 
approximation  to  the  interpolation  polynomial  based  on  these  points.  The  objective  is  to  recover 
exponentially  accurate  point  values  over  any  subintervaJ  [a,  6]  of  analyticity  of  /(x).  As  before, 
we  will  separate  the  analysis  of  the  error  into  two  parts:  truncation  error  and  regularization  error. 
Truncation  error  measures  the  difference  between  the  exact  Gegenbauer  coefficients  with  A  ~  A, 
and  those  obtained  by  using  the  collocation  point  values.  This  will  be  investigated  in  Section  3. 
The  regularization  error  measures  the  difference  between  the  Gegenbauer  expansion  using  the  first 
m  N  Gegenbauer  coefficients  with  A  ~  iV,  and  the  function  itself,  in  a  sub-interval  [a,  6],  in  which 
the  function  is  assumed  analytic.  This  error  is  estimated  in  [5]  and  we  will  si*^ /ly  quote  the  result 
in  Section  4.  The  results  are  summarized  in  Theorem  4.2.  Section  5  contains  a  numerical  example 
to  illustrate  our  results.  In  Section  2  we  shall  give  some  useful  properties  and  estimates,  including 
a  useful  approximation  result. 

Throughout  this  paper,  we  will  use  A  to  denote  a  generic  constant  or  at  most  a  polynomial 
in  the  growing  parameters,  as  will  be  indicated  in  the  text.  It  may  not  be  the  same  at  different 
locations. 

2  Preliminaries 


We  first  collect  some  useful  results  about  the  Gegenbauer  polynomials,  to  be  used  later.  We  rely 
heavily  on  the  standardization  in  Bateman  [1]. 

Definition  2.1.  The  Gegenbauer  polynomial  C’^(x),  for  A  >  0,  is  defined  by 


(1  -  [(1  -  x^r+*-i 


(2.1) 


where  G(A,ti)  is  given  by 


^  (-irr(A-H  l)r(n-|-2A) 
2’»n!r(2A)r(n  +  A -h  i) 


(2.2) 
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for  A  >  0,  by 


2”-inr(n  +  i) 

for  A  =  0  and  n  >  1,  and  by 

G(0,0)=1. 

for  A  =  0,  n  =  0.  Notice  that  by  this  standardization,  C^(x)  is  defined  by  (see  [1]): 

^o(x)  =  1, 


where  Tn(x)  are  the  Chebyshev  polynomials. 


(2.3) 


(2.4) 


(2.5) 


Under  this  definition  we  have,  for  A  >  0, 


C^(i)  = 


r(n  +  2A) 
n!r(2A) 


and 


The  Gegenbauer  polynomials  are  orthogonal  under  their  weight  function  (1  — 

-  x^)^-iC^(x)C^(x)dx  =  Sk,nhn 

where,  for  A  >  0, 


(2.6) 


(2.7) 


(2.8) 


(2.9) 


T(A)(n  +  A) 

We  will  need  to  use  heavily  the  asymptotics  of  the  Gegenbauer  polynomials  for  large  n  and  A. 
For  this  we  need  the  well-known  Stirling’s  formula: 


(27r)2x*+2c“*  <  r(x  -f  1)  <  (27r)ii"=+*  e-*eli?,  x  >  1 
Lemma  2.2.  There  exists  a  constant  A  independent  of  A  and  n  such  that 

The  proof  follows  from  (2.9)  and  the  Stirling’s  formula  (2.10). 


(2.10) 


(2.11) 


•0IVi 


□ 


Let  us  assume  that  the  function  f{x)  is  analytic  in  [-1,1]  and  satisfies  the  following 
Assumption  2.3.  There  exists  a  constant  p  >  1  and  C(p)  such  that,  for  every  A  >  0, 

max 


(x) 


< 

r 


(2.12) 

□ 


This  is  a  standard  assumption  for  analytic  functions,  p  is  the  distance  from  [-1,1]  to  the  nearest 
singularity  of  /(x)  in  the  complex  plane  (see  for  example  [7]).  For  such  function  /(x)  we  have  the 
following 


Lemma  2.4.  If  /(x)  is  an  analytic  function  satisfying  the  Assumption  2.3,  then  for  any  integer  m 
such  that  0<Tn<A— ^,we  have 

<C(rtr(A+i)  (2+i)"'.  (2,13) 


max 

-1<X<1 


dx^ 


[(1  -  x^)»-i/(x)] 


Proof; 


dx' 


(l-x^)‘-i/(x)]  = 


E 

k=0 


ml 


i—k 


^  fc!(7n  —  fc)! 


rr/(^) 


(1  -  x“)*"5 


<  cwE 


ml  (m-Jt)!2‘«r(2(A-»;))r(A+l) 


(1-x') 


2\\-k-k 


^  fc!(m  -  k)l  Y[X-k+l)  r(2A  -  k) 

m!  (m-A:)!2*'A;!r(2(A-fc))r(A+i)  r(2A  -  A:) 


^  r  (a  -  it  +  i)  r(2A  -  ifc)  A:!r  (2(A- A;)) 


2p 

m—k 


<  C(.)2™r(A.l)|;^^^^^(i 

<  C(.)2"T(A+i)|;^^(±) 

=  C(.)r(A+i)(2+ J 


m—k 


where  we  have  used  the  analyticity  assumption  (2.12),  and  the  formulas  (2.1)-(2.2),  for  the  second 
inequality.  We  have  also  used  (2.6)-(2.7),  and  the  fact  that  (l—x^)^~^~^  <  1  because  A  — A:-i  >  0, 
for  the  third  inequality.  In  the  fourth  inequality  we  have  used  the  fact  that  T  ^A  —  A*  +  >  (rn  — A)! 

because  A  —  A+|>m  —  A+l. 
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Finally,  we  would  need  to  quote  some  approximation  results  for  the  interpolation  polynomials. 

If  the  point  values  /(X|)  of  the  function  f{x)  are  known  at  the  N  Gauss  or  Gauss-Lobatto  points 

{it}  of  the  basis  functions  {^;t(a:)},  we  will  denote  the  unique  interpolation  polynomial  by 

N 

InUx)  -  XI  fkMx)  (2.14) 

it=i 

which  satisfies 

lNfixi)=^  fixi),  (2.15) 

We  then  have  the  following  error  estimates  for  the  interpolation  polynornial  I[^f(x): 

Lemma  2.5.  If  f(x)  has  m  continuous  derivatives  in  [—1,1],  then  the  interpolation  polynomial 
Jn/Ix)  defined  by  (2.14),  (2.15),  where  the  basis  functions  {(^jt(i)}  are  either  trigonometric  poly¬ 
nomials  or  the  Gegenbauer  polynomials  with  n  >  —\  (which  includes,  for  example,  the 
Chebyshev  case  /i  =  0  and  the  Legendre  case  fx  =  |),  satisfies  the  following  estimate: 

(2.16) 

where  the  weighted  norm  is  defined  by 

WfWh  =  /  t/7(x)|/(i)pdz  (2.17) 

J-i 

and  i4  is  a  constant  independent  of  N  and  m. 

Proof;  The  proof  of  this  Lemma  can  be  found  in  [8]  for  the  trigonometric  polynomial  case,  and 
in  [2]  (Theorem  6.6.1,  first  line  of  formula  (6.6.3),  for  which  the  only  restriction  about  the  weight 
corresponds  to  our  /i  >  —  5)  for  the  Gegenbauer  polynomial  case.  Notice  that  in  both  situations 
the  constant  A  is  actually  independent  of  m,  although  the  authors  did  not  explicitly  point  this  out. 

□ 

3  Truncation  Error 

Consider  an  L\  function  f{x)  defined  in  [-1,1].  We  assume  that  the  point  values  /(i,)  at  the  N 
Gauss  or  Gauss-Lobatto  points  (i.e.,  the  zeros  of  <^;v+i(®),  or  the  zeros  of  :^<i>N{x)  plus  the  two 
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boundary  points  ±1),  are  given.  Here,  the  orthogonal  basis  {<^fc(i)}  is  either  the  trigonometric 
polynomials  or  the Gegenbauer  polynomials  with  >  -i. 

are  interested  in  approximating  the  Gegenbauer  expansion  of  /(i),  with  A  ~  AT,  based  on 
a  sub-interval  [a,  6]  C  [—1,1]  of  anaJyticity.  We  start  by  introducing  the  local  variable 


Definition  3.1.  The  local  variable  ^  is  defined  by 


(  =  ^(x)  = 


X  —  S 


X  =  x(^)  -  +  S 


where 


b  —  a 


b  +  a 


Thus  when  a<i<6,  -1<^<1. 


We  assume  that  the  function  f{x)  is  analytic  in  this  sub-interval  [a,  6]  satisfying  the  Assumption 


max 

a<*<6  dx*'  ^  ^ 


for  any  A:  >  0,  with  some  p>  1. 


The  function  /(i)  has  a  Gegenbauer  expansion  in  the  sub-interval  [a,  6],  with  A  ~  A.  With 


€  and  6  defined  in  (3.1)-(3.2),  we  have 


/k  +  ^)  =  E/"(oq\o, 

/=0 

where  the  Gegenbauer  coefficients  f^{l)  are  defined  by 


•1  <  e  <  1 


/,‘(0  =  ^  J'y  -  m 


Our  goal  in  this  section  is  to  find  a  good  approximation  to  the  first  m  N  Gegenbauer 
coefficients  f^{l)  in  (3.5). 

BdiSed  on  the  known  point  values  /(xi),  we  define  not  the  usual  interpolation  polynomial  In  fix) 


in  (2.l4)-(2.15),but 


gN(^)  =  ^n(«-  /)(*)> 


6 


where  the  function  tt(x)  is  defined  by 


a  <  X  <b, 
otherwise. 


(3.7) 


with  ^(x)  defined  by  (3.1).  This  works  for  the  case  /^j  <  oo,  where  w(x)  is  the  weight 

function  of  the  basi.s  Most  commonly  used  bases  belong  to  this  class,  for  example  the  Fourier 

case  where  w{x)  =  1  and  the  Gegenbauer  case  where  tu(x)  =  (1  —  x'^)^~^ ,  when  fJ.  <  ^.  Otherwise, 
if  the  basis  {<i>k}  consists  of  the  Gegenbauer  polynomials  C^(x)  with  //  >  |,  we  should  choose 


/  (I- 


a  <  X  <  b, 
otherwise. 


Intuitively,  the  function  gf\/(^x')  has  about  X  N  continuous  derivatives,  hence  would  produce 


nice  error  estimates. 


Our  candidate  for  approximating  the  Gegenbauer  coefficients  in  (3.5)  is: 


iHi) = £^snU( + 


if  (3.7)  is  used  to  define  a(x),and  is 


g^{l)=  ^  j  -  x{^f)i~*gNi(^  + 


if  (3.8)  is  used  to  define  afx). 

How  well  do  5j^(/)  approximate  To  answer  this  question  we  define 

Definition  3.2.  The  truncation  error  is  defined  by 

m 

TE{X,  m.N,,)=  mj^  \  "  j^(0)Ci'W)l 

j_o 

where  f£(l)  are  defined  by  (3.5)  and  §^(1)  are  defined  by  (3.9). 


(3.9) 


(3.10) 


(3.11) 


The  truncation  error  is  the  measure  of  the  distance  between  the  true  Gegenbauer  expansion  in 
the  interval  [a,i»)  and  its  approximation  based  on  the  collocation  point  values  of  /(x)  in  [-1,1]. 
Using  the  result  of  the  previous  section,  we  have  the  following  estimate  for  the  truncation  error: 
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Lemma  3.3.  The  truncation  error  can  be  estimated  by 


TE{X,in,  N^()  <  A 


m(m  4-  A)r  r(Tn  +  2A)  j- 1  1  ^  j  A 


\/Am!r(2A)Ar^-|-t 

PrpoL  We  shall  just  prove  the  case  when  (3.7)  and  (3.9)  are  used.  For  this  case  we  have: 

m 

TE{X.m,N,i)  <  EI/'W-A^OIC/CI) 

1=0 

s  E  [(1  -  + «)  -  + <)]  ci'im 

^  (C/^(l))^  1  /-i 

<  12  ••  ••  ;,>  ^  -  \aix)f(x)-gN{x)\dx 

(c'j^(l)^  //■I 

£  E  (/_,  <*")  (I, 


(3.12) 


^(g^l))'  /t  „ 


/=0 


[>4] 
bn  E 


t'o 


r  A  + 


<  A 


m(m  +  A)r  (a  +  i)  r(Tn  +  2A) 


J 


E 

(=0 


(gi'(i))' 
ehf 


vOtm!  r(2A);V^-| 

where  we  have  used  (2.7)  for  the  first  inequality;  the  definitions  of  f^{l)  and  5,^(/),  in  (3.5)  and  (3.9), 
for  the  second  inequality;  the  definition  of  o;(x)  in  (3.7)  for  the  third  inequabty;  Cauchy-Schwartz 
inequality  for  the  fourth  inequabty;  the  approximation  result  (2.16)  for  the  fifth  inequabty,  with 
[A  —  “]  denoting  the  greatest  integer  <  A  —  A;  (2.13)  for  the  sixth  inequabty;  and  (2.6)  and  (2.11) 
for  the  last  inequabty. 

The  case  when  (3.8)  and  (3.10)  are  used  instead  of  (3.7)  and  (3.9)  can  be  proven  in  a  similar 
fashion.  We  only  need  to  notice  that  the  function  ^  always  analytic  in  a  <  x  <  6, 

hence  can  be  combined  with  /(x)  in  the  analysis.  This  accounts  for  the  extra  factor  N~i  in  the 
denominator. 


□ 
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We  are  now  ready  for  the  main  theorem  of  this  section: 

Theorem  3.4.  Let  the  truncation  error  be  defined  in  (3.11).  Let  A  =  ocN  and  m  =  with 
0  <  <  1,  then 


TE{iXiN,(iiN,N,€)  <  A  .  (-  +  —Y 

^  \  e"a“/3^  V2  4p/ 

where  A  grows  at  most  as  iV  2  .In  particular,  if  a  =  ;9  <  — 

TE{a(N,a(N,N,€)  < 


t;V 


where 


Proof;  The  theorem  follows  from  Lemma  3.3  and  the  Stirling’s  formula  (2.10). 


(3.13) 


(3.14) 


(3.15) 


□ 


We  now  make  two  remarks: 

Remark  3.5.  Unlike  in  the  Galerkin  case  [3],  [5]  and  [6],  where  for  the  purpose  of  truncation 
error  being  exponentially  small  we  only  need  to  assume  that  the  function  f{x)  is  in  ,  here  for 
the  collocation  case  we  need  to  assume  that  the  function  f{x)  is  actually  analytic  in  the  relevant 
sub-interval  [a,  6]. 

□ 

Remark  3.6.  We  would  like  to  point  out  that  the  collocation  procedure  (3.6),  where  we  collocate 
not  the  function  f{x)  but  its  product  with  a  weight  (1  — x^)'^~?,  is  somehow  nonstandard.  Another, 
probably  more  natural,  procedure  of  collocation,  is  to  interpolate  f{x)  with  ////(a:)  and  then  define 
the  approximate  Gegenbauer  coefficients  g^{l)  by 

9."(o  =  A  /_' 0  -  + i)chodi.  (3->6) 

instead  of  (3.9).  Unfortunately,  we  have  been  unable  to  prove  that  the  truncation  error  is  expo- 
nentijilly  small  for  (3.16).  In  fact  it  can  be  shown  that  for  some  cases  the  aliasing  error  reduces 


Iho  accuracy.  However  luuncrical  experiments  in  Section  5  seem  to  indicate  that  both  procedures 
produce  similarly  good  results. 

□ 

4  Regularization  Error  and  the  Main  Theorem 

The  second  part  of  the  error,  which  is  called  the  regularization  error  and  is  caused  by  using  a 
finite  Gegenbauer  expansion  based  on  s  sub-interval  [a,  6]  C  [— T  1]>  to  approximate  a  function  /(x) 
which  is  assumed  analytic  in  this  sab-interval,  has  been  studied  in  [5].  We  will  thus  just  quote  the 
result. 

We  again  assume  that  /(x)  is  an  analytic  function  on  [a,  6]  satisfying  the  Assumption  2.3,  see 
(3.3). 

Let  us  consider  the  Gegenbauer  partial  sum  of  the  first  m  terms  for  the  function  f{(^  +  ^): 

m 

/=0 

with  f  and  6  defined  by  (3.1)  and  (3.2),  and  the  Gegenbauer  coefficients  based  on  [a,b], 

defined  by  (3.5). 

The  regularization  error  in  the  metximum  norm  is  defined  by: 

m 

m,  0  =  /(f{  +  <)  -  E  f'WCfU)  (4.2) 

/=:0 

We  have  the  following  result  for  the  estimation  of  the  regularization  error,  when  A  ~  m  [5]: 


Theorem  4.1.  Assume  A  =  'ym  where  7  is  a  positive  constant.  If  /(x)  is  analytic  in  [a,  6]  C  [—1,1] 
satisfying  the  Assumption  2.3  in  (3.3),  then  the  regularization  error  defined  in  (4.2)  can  be  bounded 
by 


RE{jm,m,()  <  Aq^ 

(4.3) 

where  q  is  given  by 

e(l  +  27)’+27 
^  ~  /92i+2'>7'r(l  -f  7)l+'r 

(4.4) 
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which  IS  always  less  than  1.  In  particular,  if  7  =  1  and  m  =  piN  where  /3  is  a  positive 
then 


with 


RE{p(N,f3(N,e)  < 


f27(Y 


constant, 


(4.5) 


(4.6) 


□ 


We  can  now  combine  the  estimates  for  truncation  errors  and  regularization  errors  to  obtain  the 
following  main  theorem  of  this  paper: 

Theorem  4.2.  (Removal  of  the  Gibbs  Phenomenon  for  the  collocation  case). 

Consider  a  Lj  function  f(x)  on  [—1,1],  which  is  analytic  in  a  sub-interval  [a, 6]  C  [—1,1]  aud 
satisfies  the  Assumption  2.3  in  (3.3).  Assume  that  the  point  values  /(a;,)  at  the  N  Gauss  or  Gauss- 
Lobatto  points  of  {<^fc(a:)},  are  given.  Here  {<f>k(x)}  is  either  the  trigonometric  polynomials  or 
the  Gegenbauer  polynomials  Cj^(x)  for  /j.  >  — Let  0  <  /  <  m  be  the  Gegenbauer  expansion 
coefficients,  defined  in  (3.9)  or  (3.10),  based  on  the  sub-interval  [a,  6],  of  the  collocation  polynomial 
gN{x)  defined  in  (3.6).  Then  for  A  =  m  =  P(N  with  /3  < 


where 


max 


qr  = 

rLta 


rn 


1=0 


<  A  (^q^  +  q)i  ^ 


(4.7) 


27(3 

f>  +  -)] 

.  2e 

V  2pJi 

<  1, 


QR 


f27(Y  , 

V32^j  ^ 


and  A  grows  at  most  os  N  2 
Proof:  Just  conbine  the  results  of  Theorems  3.4  and  4.1. 


We  remark  that  no  attempt  has  been  made  to  optimize  the  parameters. 


ll 


5  Numerical  Results 


In  this  section  we  give  a  numerical  example  to  illustrate  our  result.  We  will  test  both  the  Fourier 

and  the  Chebyshev  cases. 

Example  5.1.  We  take  the  following  function 

ffx)  =  /  sin(cos{a;)),  if  a<x<b 
^  1  0,  otherwise 

and  try  to  recover  pointwise  values  of  f{x)  for  a  <  x  <  6.  This  function  is  also  used  in  [6].  In  order 
to  see  the  effect  of  relative  locations  of  the  sub-interval,  for  the  Chebyshev  case  we  will  consider 


two  situations  [a,  6]  =  [-0.5, 0.5]  and  [a,  6]  =  [0,1].  For  the  Fourier  case  we  will  only  consider  the 
case  [a,  b]  =  [-0.5, 0.5]  because  of  periodicity. 

We  first  consider  the  Fourier  case.  We  assume  that  the  point  values  /(x.),  for  the  21V -I- 1  equaUy 


spaced  points: 


2i 


Xi 


2N  4-  1 


-N  <i<  N, 


(5.2) 


are 


given.  Two  coUocation  procedures  are  tested.  Procedure  A,  which  is  the  procedure  proven  in 
Section  3  to  be  exponentially  accurate,  is  to  interpolate  a{x)f{x)  as  in  (3.6)  and  (3.7),  and  then 
build  the  approximating  Gegenbauer  coefficients  g^{l)  by  (3.9),  finally  use 


in 

1=0 

as  our  approximation  to  /(x)  in  [a,  6].  We  also  test  the  more  natural  collocation  Procedure  B, 
for  which  we  interpolate  /(x)  as  in  (2.14)-(2.15),  and  then  build  the  approximating  Gegenbauer 
coefficients  g^il)  by  using  (3.16),  finally  using  (5.3)  to  approximate  /(x)  in  [a,  6]. 

For  the  Fourier  method,  since  the  collocation  points  are  equally  spaced,  we  also  have  the  choice 
to  use  all  the  coUocation  points  in  [-1,1]  to  build  the  approximation  in  [a,  6]  (referred  to  as  global- 
local)]  or  use  just  the  coUocation  points  inside  [a,  6],  with  suitable  scaling  to  expand  [a,  6]  to  [-1, 1] 

(referred  to  as  local-local). 

We  shaU  use  the  parameters 


m  =  0.2(N,  A  =  0.4ciV 


(5.4) 
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where  2N  +  1  is  the  total  number  of  collocation  points  in  the  full  interval  [-1, 1].  This  choice  is  also 
used  in  [6]  (with  the  right  scaling).  We  will  use  (5.4)  for  all  the  test  cases  of  Fourier  collocation. 
No  attempt  has  been  made  to  optimize  these  parameters  for  each  individual  case. 

In  Fig.  1  we  show  the  errors,  in  a  logarithm  scale,  for  =  10,  20,  40  and  80,  for  the  Procedure 
A  (global'local  at  left,  local-local  at  right).  In  Fig.  2  we  show  similar  errors  for  the  Procedure 
B.  We  can  clearly  see  that  both  Procedure  A  and  Procedure  B  produce  similarly  good  results. 
Also,  it  seems  that  global-local  and  local-local  are  producing  similar  errors  for  the  same  choice  of 
A  and  m.  Considering  that  the  effective  N  (the  number  of  collocation  points)  is  smaller  for  the 
local-local  case,  hence  round  off  errors  should  also  be  smaller,  it  seems  that  local-local  should  be 
used  if  possible  (see  also  the  discussions  in  [9]).  However,  there  are  situations  (for  example  when 
we  are  solving  a  PDE  and  is  obtained  from  the  time  evolution)  where  the  point  values  /(x,) 
may  not  be  accurate,  however  the  interpolation  polynomial  In/  (which  may  interpolates  f{x)  at 
some  other  points  than  x,)  is  accurate.  In  such  situations  global-local  is  more  appropriate. 


Fourier,  Global-Local,  Collocate  f*welght 


Fig.  1:  Errors  in  log  scale,  /(x)  defined  by  (5.1).  [a,  6]  =  [—0.5, 0.5].  Fourier  collocation  with 
Procedure  A,  namely  interpolating  /(x)  *  a(x)  with  q(x)  defined  by  (3.7).  A  =  0.4(W  and  m  = 
0.2(N,  N  =  10,20,40,80.  Left:  global-local;  Right:  local-local. 
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Fig.  2:  Errors  in  log  scale,  f(x)  defined  by  (5.1).  [a,  6]  =  [—0.5, 0.5].  Fourier  collocation  with 
Procedure  J9,  namely  interpolating  /(x)  only.  A  =  OAeN  and  m  =  0.2€iV,  N  =  10,20,40,80.  Left: 
global- local]  Right:  local-local. 

Next  we  consider  the  Chebyshev  case.  We  assume  that  the  point  values  /(li),  at  the  +  1 
Gauss- Lobatto  points 

X,-  =  cos  >  0  <  i  <  N,  (5.5) 

are  given.  We  again  test  both  Procedure  A  and  Procedure  B,  now  for  the  middle  interval  [a, 6]  = 
[—0.5, 0.5]  as  well  as  for  the  boundary  interval  [a,  6]  =  [0, 1].  The  choice  of  A  and  m  are  the  same 
as  in  [6]: 

m  =  O.lfAT,  A  =  0.2(N  (5.6) 

This  is  also  the  same  as  (5.4)  for  the  Fourier  case,  considering  that  now  the  number  of  collocation 
points  IS  N  +  I  instead  of  2N  -f  1. 

In  Fig.  3  we  show  the  errors,  in  a  logarithm  scale,  for  N  =  20,  40,  80  and  160,  for  the 
Procedure  A  ([a,  6]  =  [—0.5, 0.5]  at  left,  [a,  6]  =  [0,1]  at  right).  In  Fig.  4  we  show  similar  errors 
for  the  Procedure  B.  We  can  again  see  that  both  procedure  A  and  Procedure  B  produce  similarly 
good  results.  Also,  these  errors  are  comparable  with  those  obtained  in  [6]  where  Galerkin  instead 
of  collocation  is  used. 
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Fig.  3:  Errors  in  log  scale,  f{x)  defined  by  (5.1).  Chebyshev  collocation  with  Procedure  A,  namely 
interpolating  f{x)*a{x)  with  a(x)  defined  by  (3.7).  A  =  0.2€N  and  m  =  O.UN,  N  =  20,40,80, 160. 


Left:  [a,  6]  =  [-0.5, 0.5];  Right:  [d,b]  =  [0, 1]. 


Fig.  4:  Errors  in  log  scale,  /(i)  defined  by  (5.1).  Chebyshev  collocation  with  Procedure  B,  namely 
interpolating  f{x)  only.  A  =  0.2(N  and  m  =  0.\(N,  N  =  20,40,80,160.  Left:  [a, 6]  =  [—0.5, 0.5); 
Right:  [a,b]  =  [0,1]. 

6  Concluding  Remarks 

We  have  proven  the  exponential  convergence  in  the  maximum  norm,  of  a  reconstruction  procedure 
using  Gegenbauer  series  based  on  C'j^(x)  with  large  A,  for  any  L\  function  in  any  sub-interval  [a,  6] 
in  which  the  function  is  analytic,  if  we  are  given  the  collocation  point  values  of  the  function  at  the 


Cia\iss  or  Gauss-Lobatto  ])oints.  The  proof  covers  Fourier,  Chebyshev,  Legendre  and  more  general 
Gegenbauer  collocation  methods.  Numerical  examples  are  also  given. 
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